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I Abstract. We study average growth of the spectral function of the Lapla- 

cian on a Riemannian manifold. Two types of averaging are considered: with 
respect to the spectral parameter and with respect to a point on a manifold, 
^-li We obtain as well related estimates of the growth of the pointwise ^-function 

along vertical lines in the complex plane. Some examples and open problems 
regarding almost periodic properties of the spectral function are also discussed. 
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^ ■ 1. Introduction 

- 

■ 1.1. Notation. Throughout the paper, we use the following notation. 

• C denotes various positive constants whose precise values are not impor- 
tant for our purposes; dependence on various parameters is indicated by 

^ ! lower indices. 

I • dx^y is the Riemannian distance between x and y. 

• M+ := (0,+oo). 

• t,s,X,fi,TeR and ^ G M'^. 

• A+ := max{A,0}, [A] = min{n G Z| A ^ n} is the ceiling function and 
2 I [AJ = max{n G Z| A ^ is the floor function. 

(yQ I • /(^) = J /(x)e^*^^ dx is the Fourier transform of / and (Z)"^ is the inverse 

^ ' Fourier transform, so that (/)^ = /. 

• * denotes the convolution on M. 

^ . • r is the gamma-function and [^) are the binomial coefficients. 

^ i • Jo is the Bessel function of the first kind of order a. 

Let / and g be real- valued functions on M+, and let > 0. We write 

. /(A) = Oi\P) if limsup,^+^ A-P/(A) < oo; 
. /(A) = o{Xn if lim,_.+„, A-P/(A) = 0; 

• /(A) = 0(A-°°) if A"/(A) ^ as A ^ +00 for all m G M+ ; 

• o{g{X)) if there exists a sequence /x„ +oo such that ^ 
C > for all n; 

• /(A) 2> g{X) if there exists Aq such that /(A) > C g{X) for all A > Aq. 
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1.2. Spectral function. Let M be a compact Riemannian n-dimensional man- 
ifold and A be the Laplacian on M with the eigenvalues = Aq < ^ A2 ^ 
. . . Aj ^ ... and the corresponding orthonormal basis of eigenfunctions {(f) j{x)}. 
Let Nx,y{X) be the spectral function, 



:i.2.1) iV.,,(A) := 



0<Aj<A 

The spectral function N^^y^X) is the integral kernel of the spectral projection of 
the operator \/A, corresponding to the interval (0, A) . This way the spectral 
function could be defined in the non-compact case as well, see Example 11.3.11 
By the Weyl formulae, 

(1.2.2) N^JX) = s + 0(A"-^) 

and 

(1.2.3) iV,,,(A) = 0(A"-i) 

for all a; ^ y G M (see [SI Ei lEl]. The estimates (11X21) and ([L23D are 
sharp and attained on a round sphere. However, the bounds can be improved to 
o(A"'~^) instead of 0(A""^), provided that the set of geodesic loops originating at 
X and the set of geodesies joining x and y are of measure zero (see [Si] or jSVj ). 

The following condition is used to prove the lower bounds on the spectral 
function (see subsection 15. ip . 

Condition 1.2.4. The points x and y are not conjugate along any shortest 
geodesic segment joining them. 



It follows immediately from [Mill Corollary 18.2] that for each x G M, Con- 



dition [L231 is true for almost all G M. As was shown in |JP[ Theorem 1.1.3] 
that Condition 11.2.41 implies 



:i.2.5) N,,y{\) ^ o( 



A^ 



1.3. Discussion. Let cF^^yit) be defined by {ax,y{t)y = (j)j{x) 4>j{y) 6{X — Xj) . 
By the spectral theorem, (Tx,y{t) coincides with the distributional kernel of the 
operator exp(— zt\/A) or, in other words, with the fundamental solution of the 
corresponding hyperbolic equation. Proofs of fll.2.2l) - fll.2.5l) are based on the 
study of singularities of the distribution ax,y(t), which is usually done by means 
of Fourier integral operators (see, for example, [DC] or |SV] ). It is well known that 
these singularities lie in the set Cx,y of lengths of all geodesic segments x,y E M. 
In particular, CTx,y{t) with y ^ x is infinitely smooth in a neighbourhood of the 
origin, whereas crx,x{t) has a strong singularity at t = related to the main term 
in the Weyl formula fll.2.2p . If Condition 11.2.41 is fulfilled then the singularity at 
t = dx^y can be explicitly described, which leads to the lower bound fll.2.5p . 
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The purpose of this paper is to study the average growth of Nx^y{X), where 
averaging is considered either with respect to the spectral parameter A or with 
respect to a point y. The following examples indicate that, while both estimates 
(11 .2. 31) and (11. 2. 51) are sharp, the typical order of growth of the spectral function 
is more likely to be A^ rather than A*^"^. 

Example 1.3.1. (cf. |Pet Example 4.1]) For the Laplacian in the spectral 
function is given by 

1 r \n 

where xb" is the characteristic function of the unit ball. Taking into account 
[Pel formula (4.4)] and [GRl formula 3.771(8)], one gets N:,^y{\) = U{X,d^,y) 
where 

(1.3.2) Ar(A,4J := (27r)-trfx;|At J|(4,,A) 

and dx,y = \x — y \ . Using [GR[ formula 8.451(1)] we obtain for all x ^ y: 

Example 1.3.4. Let S" be a round sphere of dimension n ^ 2. Given x G 
denote by —x the diametrically opposite point. Then Nx.y{X) = O ^A^ j for all 

y^±x and iV.,_.(A) ^ o(A"-i). 

1.4. Plan of the paper. The paper is organized as follows. In the next section 
we present our main results. In Section [3] we discuss some examples and open 
problems regarding almost periodic properties of the spectral function. In Section 
m we give proofs of the upper bounds formulated in the subsections 2.1-2.4. 
Proofs of the lower bounds stated in the subsection 2.5 are given in Section [51 
All examples are justified in Section [61 Finally, in Section [7| we establish the 
estimates on the (^-function presented in the subsection 2.6. 

2. Main results 

2.1. Average over the manifold. The principal object of study in the present 
paper is the rescaled spectral function 

~ 1 — ji 

(2.1.1) Nx^y{\) := A— Ar,,^(A) 

on a compact n-dimensional Riemannian manifold M. Note that, by (11.3.31) . in 
the Euclidean case the rescaled spectral function J\f{X,dx^y) := Af{X, d^^y) 
is bounded for each fixed x ^ y . Our main result is 
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Theorem 2.1.2. For any compact n- dimensional Riemannian manifold M there 
exists a constant Cm such that 



2 

^ Cm, Vx g M, VA e M+, 



(2.1.3) / N,^y{\)-N{X,d,^y 

where d^ y is the Riemannian distance between x,y E M . 



Remark 2.1.4. Note that Nx,x{X) ^ A^. We subtract the term Af{\,dx,y) in 
(12.1.31) in order to "regularize" the rescaled spectral function in a neighbourhood 
of the diagonal. Geometrically this can be interpreted as follows. Consider a 
normal coordinate system centred at x. In these coordinates the Riemannian 
metric is Euclidean at the point x, and hence near x the geometric Laplacian 
A can be viewed as a perturbation of the Euclidean Laplacian (a similar idea 
was used in |Polt Section 3] and goes back to [AK' Theorem 6.1]). We justify 
the regularization using the Hadamard parametrix for the wave kernel, see the 
subsection 14.31 

2.2. Spectral average. Theorem 12.1.21 implies a number of results. The first 
one describes the growth of the spectral function on average with respect to A. 

Theorem 2.2.1. For every finite measure v on ]R_|_ and each fixed x G M, there 
exists a subset M^^u C M of full measure such that 



/■oo 

(2.2.2) / |iV,,j,(A)|Mz/(A) < oo, Vy G 

Jo 



It should be emphasized that, generally speaking, the set M^^i, depends on the 
choice of z/. Otherwise fl2.2.2p would mean that Nx,yW is bounded for almost 
all y E M , which is not always the case (see the next subsection). 

In particular, taking du = (A + 1)~^ (In A)~^~^ dA and applying Theorem 12.2. 1^ 
we see that 



(2.2.3) / A-^(lnA)-^~"|iV^,j,(A)pdA < oo 

Jo 

for all e > 0, all X G M and almost all y G M. 



Remark 2.2 A. This paper was inspired by [Ran], where the estimate fl2.2.2p was 
proved for surfaces of constant negative curvature under the assumption that 
dz/(A) = (A + 1)^^^^ dA with some e > 0. Randol's proof is based on the estimate 
(12.6.41) for the pointwise ^-function of the Laplacian and uses some results of 
complex analysis obtained in |HIP] . We use a more direct approach that is 
applicable in higher generality and gives better estimates. It also allows us to 
improve the bounds for the (^-function. Even though these estimates are not 
needed in our proof, they seem to be of independent interest. We have included 
them in Section 12. 6^ where we also explain the relation between Theorem 12.2.11 
and properties of the ^-function. 
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2.3. Growth of the rescaled spectral function. Theorem 12.2.11 and Exam- 
ple [L3]ll suggest that Nx^y{X) = 0(1) for almost all x,y E M on any manifold 
M. However, this is not true. In particular, for any negatively curved manifold 
M there exists a constant um > depending on certain dynamical properties of 
the geodesic flow, such that Nx,yW 7^ o ((In for all x,y e M (see [JP]l At 
the same time, for flxed x and y on a negatively curved manifold, the sequence 
of A-s yielding the logarithmic growth of the rescaled spectral function is very 
scarce. This sequence is quite sensitive to the choice of the points x and y. In 
particular, for flxed x and A, the set of points y, for which the function Nx,y{X) 
has a logarithmic "peak" , is expected to be very small. 

The theorem below shows that a similar effect takes place on any manifold, 
not necessarily negatively curved. 

Theorem 2.3.1. There exists a constant Co not depending on M , and a constant 
Cm depending on the geometry of M, such that for any C > Cq, the measure of 
the set 

n + l 

(2.3.2) ^^.(A,^) := {yeM : |iV,,j,(A)| ^ C (/i + rf.,, ^ )} 

satisfies 

meas(f22;(A, /i)) ^ *^ 



C^fx^ 

for any point x E M and any A G M+. 

Corollary 2.3.3. Let {^k} CLnd {r^} he positive increasing sequences converging 
to +00. Then, for each fixed x G M , the measure of the set of points y E M such 
that 

(2.3.4) \N.,y{n)\ > fik, Vfc = 1,2,..., 
is equal to zero. Moreover, if /i^^ < oo then 

(2.3.5) limsnp\fi^'N,,y{Tk)\ = 

for almost all y G M. 

2.4. Weighted average over the manifold. Let us consider the average of the 
rescaled spectral function over M with a weight given by a power of the distance 
function d^^y 

Theorem 2.4.1. For all h ^ 0, x E M and A G IR+, we have 

C^,m(1 + A1-'^), x^l. 



(2.4.2) / rf-JiV.,,(A)|2dy ^ 



CAf (l + |lnA|) 



X 



Note that J^^ |iV^,j/(A)P dy = A^"" N^^^iX) = CX + 0(1) . Therefore the esti- 
mate (12.4.21) with X = is order sharp. We also remark that, as follows from the 
proof of Theorem 12.4.11 the constant C^^m blows up as pj-^ when x — > 1. 
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Corollary 2.4.3. Let u{x,y) be a function on M x M such that \u{x,y)\ ^ 
C d^y with some nonnegative constants h and C , and let Kx be the operator 

defined by the integral kernel lC\{x,y) := u{x,y)\Nx^y{\)\'^ ■ Then Kx maps 
LP{M) into LP{M) and 

I^Ca/ (1 + I lnA|) , x=l. 

In particular, Corollary 12 . 4. 31 implies that the commutator of the operator given 
by the integral kernel |iV^^j,(A)p with the multiplication by a smooth function / 
is bounded in L^(M) and its norm is O (| In A|) . Indeed, the integral kernel of this 
commutator coincides with (/(x) — f{y)) \Nx,y{X)\'^ , and \f{x) — f{y)\ ^ C dx,y 
with some C > 0. 

2.5. Lower bounds. The following theorem shows that our upper estimates 
cannot be significantly improved. 

Theorem 2.5.1. // Condition \\.2A\ is fulfilled then 

(2.5.2) A-^-i / /i«|iV,,^(/i)|^d/x > 1, Vg^O, p^l. 

Jo 

In particular, A~^ |iV^._y(/i)| d/i S> 1 . 

Corollary 2.5.3. Let Condition \\.'lA\ be fulfilled, and let f be a positive function 
on M^^ . // there exists a constant g ^ such that 

(2.5.4) limsup (/i^+Mnf (r-VM) ) >0 

then J f{fi) \Nx,y{tJ')\^ dfi = oo for all p ^ 1. In particular, we have 

/■oo 

(2.5.5) / /i-i|iV,,j,(/i)|^'d/i = oo, Vp^l, 
Jo 

and 

(2.5.6) |iV,,,(/i)rd/i = oo, Vp^l, 



for every increasing sequence fik +oo. 

As follows from (|2X5|) with p = 2, Theorem EXH fails for dz/(A) = (A + 1)-MA 
on any manifold M. 

Remark 2.5.7. Theorem 12.5.11 improves upon |JPl Theorem 1.1.3]. It is quite 
possible that Condition 11.2.41 in this theorem can be removed: our proof works 
whenever cr^,y(t) has a sufficiently strong singularity. It is hard to imagine the 
situation where this does not happen; we are not aware of any counterexamples. 
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2.6. The (^-function. The function Z^y of complex variable z = t + is defined 
by 

noo POO 

(2.6.1) Z,,y{z) := / X-'dN,,yiX) = z / A"^"^ iV,,j,(A) dA 

Jo Jo 

is said to be the pointwise (-function of the Laplacian. It is the integral kernel 
of the pseudodifferential operator A" 2. It is well known that Zx^y{z) is an entire 
function on C for all fixed x ^ y (see, for example, [Sht Theorem 12.1]). 

Remark 2.6.2. Further on we call Z^^y(z) simply the (-function. Note that Z^^y{z) 
should not be confused with the function Z(z) = Tr A^i that is usually referred 
to as the (^-function of the Laplacian. 

The second equality in (12.6. ip implies that Z^^y^z) = z (J^ Nx^y){—z) where 
M is the Mellin transform, {Mf){z) := X^~^ f (X) dX . Recall that, for each 
if: G M, the Mellin transform {A4f){t + is) of a distribution / on M+ coincides 
with the inverse Fourier transform of the distribution e*'*/(e'^) on M modulo the 
factor (27r)^^. The inversion formula reads 

f{X) = {2m)-' / X-\Mf){z)dz, 

J t—ioo 

where the integral is understood in the sense of distributions. 

Obviously, e^^'f{e^') G ^^(M) if and only if A*-5/(A) G L2(M+) . Therefore 
{t + is)-^ Z^^y{t + is) G L2(R) if and only if X-^-^N.^^y{\) G L'^{^+) for each 
fixed t, and 

rt+ioo 

(2.6.3) N^^yiX) = {2m)-' X' z-' Z,^y{z) dz 

J t—ioo 

in the sense of distributions. In particular, if (to + is)~' Zx:^y(tQ + is) G L^(]R) for 
some to e R then {t + is)''' Z^^y{t + is) G L'^{R) for all t > to • 

In [Ran], for a surface of constant negative curvature, it was shown that 

Tl — 1 

(2.6.4) {t + f)-'Z,^yit + i-)eL\R), Vt>^— , 

almost everywhere. By the above, this inclusion is equivalent to Theorem 12.2.11 
with du= {X + l)-'-'dX. 

Let (s) := (1 + |sp)^/^. In the last section we shall prove the following two 
theorems. 

{Ct, if n <t, 

a(kr-*+rf*7(s)t-*) , t<f . 
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Theorem 2.6.6. For all x G M and all e > , we have 

In view of Fubini's theorem, fl2.6.7p immediately imphes fl2.6.4p for aU x and 
almost all y. However, even the improved estimate (12.6.71) does not seem to be 
sufficient to obtain our upper bounds for the spectral function. 

Remark 2.6.8. It is quite possible that Theorems 12.6.51 and 12 . 6 . 61 remain valid for 
t = n and t = f , but our proof does not work in these cases. 

2.7. Possible generalizations. All the above results can easily be extended to 
an elliptic self-adjoint pseudodifferential operator A acting on a compact manifold 
without boundary. For such an operator, one has to consider trajectories of the 
Hamiltonian flow generated by its principal symbol instead of geodesies, and to 
use Fourier integral operators or the global parametrix constructed in [SVj instead 
of the Hadamard representation (14.3.11) for the study of singularities of (Jx,yit) ■ 

Note that some of our estimates do not require the spectrum to be discrete. 
In particular, it may well be possible to extend the results that do not involve 
integration over M to the case of noncompact manifolds. For instance, the func- 
tions / pi(A — fi) dNx,y{[j) and / |pi(A — /i)p dA^s,y(yu) are the integral kernels 
of the operators pi(A — A) and |pi(A — v4)p, where A is the restriction of 
-\/A to the subspace spanned by the eigenfunctions 0i,02,---- Therefore the 
key equality (14.2. ip is easily obtained by rewriting the obvious operator identity 
(pi(A — A))* pi{\ — A) = |pi(A — A)p in terms of integral kernels. 

It would be also interesting to extend our results to an elliptic self-adjoint 
differential operator on a manifold with boundary, subject to suitable boundary 
conditions. In this case the role of geodesies is played by Hamiltonian billiards. 
One has to consider interior points x and y and to make appropriate assumptions 
to avoid problems with the so-called grazing and dead-end trajectories (see |SVj ). 



3. Almost periodic properties of the spectral function 

3.1. Besicovitch almost periodic functions. Let p ^ 1. Recall that, for a 
measurable function / on W^, its Besicovitch seminorm is defined by 

(3.1.1) WfWsP := limsupf;^ ^fWl^dX 

T^oo \^ Jo 

The space of Besicovitch almost periodic functions is defined as the completion 
of the linear space of all finite trigonometric sums Ylk=i '^fcC*^''^ with G C and 
0jt G R with respect to the Besicovitch seminorm. Clearly, B^^ C B^^ and 
II/IIbm ^ for all pi> P2- 
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For each real- valued function f E , there exists a sequence of real numbers 
dk called the frequencies of /, such that 



(3.1.2) lim \\f -y2akSmiekX + (j)k)\\Bp = 0, 

k=l 

where the coefficients G M and the phase shifts 0^ G M are some constants 
(see, for example, [Bes] ) . If f l3.1.2p holds, we shall write f sm{9kX + (j)k) ■ 

Theorem 12.2.11 motivates the following 

Conjecture 3.1.3. On any compact Riemannian manifold M , \\N.j.^y{X)\\j32 < oo 
for each fixed x E M and almost all y E M. 

Conjecture 13.1.31 holds for round spheres and flat 2-tori. In fact, in Examples 
13.2.41 and 13.2.71 we prove more than finiteness of the Besicovitch seminorm: it 
is shown that the rescaled spectral function on these manifolds is 5^-almost 
periodic for each fixed x and almost all y, which implies the finiteness of the 
Besicovitch seminorm. In both cases, the set of frequencies coincides with the set 
Cx,y of lengths of all geodesic segments joining x and y. 

Remark 3.1.4. It was proved in [BIT] and |KMSj that the rescaled error term 
in Weyl's law has an almost periodic expansion in on surfaces of revolution 
and in B^ on Liouville tori. On Zoll manifolds, the rescaled Weyl remainder 
(albeit with a different order of rescaling) has an almost periodic expansion in 
B^ |Sch] . The frequencies of these expansions are the lengths of closed geodesies. 
The spectral function, similarly to the Weyl remainder, has oscillatory behaviour, 
and hence it is natural to study it in the context of almost periodic functions. 
Moreover, as indicated by Theorem 12.2. H the order of rescaling in this case could 
be chosen universally for all manifolds of a given dimension. As was mentioned in 
Section [T73| the lengths of geodesic segments joining x and y are the singularities 
of the distribution Cx^yit), and hence play the same role for the spectral function 
as the lengths of closed geodesies for the Weyl remainder. The link between the 
spectral function and the set of lengths Cx,y has a natural interpretation from 
the viewpoint of the quantum-classical correspondence. 

3.2. Spectral function on spheres and tori. Let us start with a toy example 
— the spectral function on the unit circle S^: 

(3.2.1) iv.,(A) = iv.,(A) = i y: ^o<^^^j=-^+^ ""^^^^!;.?!'"^^ 

l^n<A 2 ^ 

Note that in dimension one no rescaling occurs, and the constant term (27r)^^ 
is subtracted, because the eigenfunction corresponding to the zero eigenvalue is 
excluded in the definition fll.2.ip . In higher dimensions the contribution of the 
constant term to A^x,y(-^) is negligible in B^ due to the rescaling. 
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Example 3.2.2. If M is the unit circle then the spectral function N^^yiX) is B^- 
almost periodic for al\x ^ y with the set of frequencies C^^y = {\dx,y + 27Tk\ }k& , 
and 

^ +00 ^ 

(3,2,3) N^JX) ~ + ^ _^^__^sm(AM„ + 2rt|). 

The next two examples generalize Example 13.2.21 to round spheres and fiat two 
dimensional tori. 

Example 3.2.4. On a fiat square 2-torus = M^/(27rZ)^, the rescaled spectral 
function is i?^-almost periodic for all x 7^ y and 

~ ^ 2 sin(A|x-?/ + 27r77| -f) 

(3.2.5) iV.,,(A) ~ (27r)3/2|x-y + 2vrr^|3/2 " 

Here | ■ | denotes the length of a vector in R^. A similar formula holds on fiat tori 
corresponding to arbitrary lattices. 

Remark 3.2.6. In dimensions n ^ 3 the situation is significantly more compli- 
cated. For example, following the argument of |Pet[ Corollary 1.4] one could 
show that ||iV'a;^j^(A)||/32 = 00 on a fiat n-torus = R"'/(27rZ)"' if the vector 
7r^^(x — y) is rational. This happens due to unbounded multiplicities in the set 
Cx^y However, it does not contradict Conjecture 13.1.31 because for each fixed 
X e T", the vector tt~^{x — y) has rationally independent coordinates for almost 
all y e T". 

Recall that the Morse index of a geodesic segment joining x and y is the number 
of points on the segment that are conjugate to x, counted with multiplicities (see 
|Milt Theorem 15.1]). Let H{x) be the "reversed" Heaviside function: H{x) = 
if X ^ and H{x) = 1 if x < 0. 

Example 3.2.7. Let M be the unit round sphere of dimension n ^ 2, and 
let X 7^ 1/ e be any two non-opposite points. Then Nx,y{X) G B"^, and 

(n— l+2a;fe)7r 



+°o 2 sin ( A y + 27rA;| 

(3.2.8) N,,y{X) ~ 



(2tc) "2 ^ (sin dx^y) "2 ^ {d^^y + 2Tik\ 

where ujk = {n — 1) {2\k\ — H{k)) is the Morse index of the geodesic segment of 
length \dx^y + 27ck\ . 



Note that, unlilce (13.2.50 . the phase shifts in the expansion (13.2.80 depend on 
the number of conjugate points on the geodesic segments. 
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3.3. An open problem. In this section we describe a possible route for gener- 
alizing Examples 13.2.41 and 13.2.71 Let M be an arbitrary compact ra-dimensional 
Riemannian manifold. Let T^^y be the set of all geodesic segments joining x and 
y. For every 7 G F^^^y, let /(7) be its length and uj{pf) be its Morse index. As be- 
fore, set Cx^y = {/(7)| 7 G F^.^^}. Along each geodesic segment 7 G Tx,y, consider 
the matrix Jacobi equation A" + RA = 0, where the coefficient R is defined in 
terms of the Riemann curvature tensor and the parallel transport along 7 (see 
|Chlt p. 104]). We assume that 7 is naturally parametrized, and A satisfies the 
initial conditions A{0,^^) = 0, A'(0,^^) = 1, where ^-y is the unit tangent vector 
to 7 at the point x. Set 0(7) = | det ^(^(7), ^^)|. For example, 0(7) = /(7)"~^ on 
a flat square n-torus and 0(7) = | sin/(7)|"~^ on a round n-sphere. In dimension 
two, 0(7) = |J(/(7))|, where J(t) is the orthogonal Jacobi field along 7 with the 
initial conditions J(0) = and J'(0) = 1. 

Problem 3.3.1. Let M be a compact n-dimensional Riemannian manifold. Is 
it true that for all x & M and almost all y & M the rescaled spectral function 
Nx,y{X) has an almost periodic expansion 

~ . . 2 V- sin(A/(7)-^^-cu(7)5) 

3.3.2 Nxy(X) ^ V — ^= 

^ ' ' (2vr)^ ^^^^ /(7)v^ 

in for some p^ 1? 

Let us show that the expansion ( ]3.3.2p is well-defined if the points x,y E M 
are not conjugate along any geodesic joining them (by (Mill Corollary 18.2], this 
condition is satisfied for any fixed x G M and almost all y G M). Indeed, if 
the points x,y are not conjugate along any geodesic, 0(7) 7^ for any 7, the set 
of lengths Cx,y is discrete, and each element has finite multiplicity, i.e. appears 
in Cx,y at most a finite number of times [Mill Theorem 16.3]. This implies, in 
particular, that the set Cx,y is infinite, because for any two points on a compact 
manifold there exists an infinite number of geodesic segments joining them |Ser j . 
If M has no conjugate points, it is easy to check that fl3.3.2p agrees with |JPt 
formula (5.1.3)]. 

Remark 3.3.3. If (13.3.21) does hold for some p ^ 1 on a given manifold, it would 
be interesting to determine the maximal possible value of p . For instance, it 
is quite likely that for round spheres one can take any p ^ 1. Note that if the 
almost periodic expansion is valid for some p > 1, then the Fourier coefficients 
of (I3.3.2P lie in Iq for q = max(2, ( |ABIl section 4]). This can be viewed as 
a dynamical condition on the manifold (see |Patl Section 3.1] for some related 
results), and it is not clear whether it always holds. At the same time, even for 
p = 1, a positive answer to Problem 13.3.11 provides a lot of useful information 
about the spectral function; in particular, it implies that the rescaled spectral 
function has a limit distribution (see |KMSt Appendix II]). Understanding the 
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properties of this distribution on a given manifold is a problem of independent 
interest . 

4. Proofs of the upper bounds 

4.1. Auxiliary functions. Let us fix a real- valued even rapidly decreasing func- 
tion p G C°°(M) satisfying the following condition. 

Condition 4.1.1. p G C^(M) , p = 1 in a neighbourhood of the origin and 

supp p C (—5, 5) for some 5 > . 

Condition 14.1.11 implies that p{t) dr = | and Jq°° p{t) dr = for all 
k = 1,2,... Let Pi(r) := sign r Jj^ p(/i) d/i if r 7^ , and Pi(0) := — | . The 
function pi is rapidly decreasing, odd and infinitely differentiable outside the 
origin r = 0, so that ^ Pi(t) = — p(t) for all r 7^ . It has a jump at the origin; 
in view of Condition HXH pi(0) = pi(-O) = -pi(+0) = -i . 

Let 

(4.1.2) N,,y.o ■■= P * N,^y , N,,y;l ■= N,^y - N,,y,o ■ 

The following elementary lemma is a slight variation of [S2i Lemma 1.2]. 

Lemma 4.1.3. N^,y-i{X) = J pi(A - fi) dN^^y{fi) for all A G M. 

Proof. If A is not an eigenvalue then, integrating by parts, we immediately obtain 

/A /»oo 
pi{X - fi) dN^^yifi) + / pi{X - fi) dN^^yifi) 
-00 J A 

= (pi(+0) - Pi(-O)) N,^y{\) - j p{\- fi) N.^yifi) d/i , 

where the right hand side coincides with Nrc^y-i{X). If A = Xj then the same 
equality holds for the function 

Ni^X) := N,,yiX) - (iV.,,(A, + 0) - N.,,yiX, - 0)) x,(A) , 

where Xj{X) is the characteristic function of the interval {Xj, 00) . Since 

poo 

Pi{Xj - fx) dxjifJ,) = pi(0) = pi(-O) = - / p(r)dr 

^0 

P(Aj - /i) Xi(/^) d/i = x{Xj) - J p{Xj - /i) Xj(/i) d/i , 
the lemma remains valid when A is an eigenvalue. □ 

Remark 4.1.4. The above lemma turns out to be very useful for obtaining esti- 
mates of the spectral and counting functions. Usually, the singularities of cT^,y(t) 
for small values of t can be described explicitly. Then, taking the inverse Fourier 
transform, one obtains full asymptotic expansion of N^^y-olX) . The asymptotic 
behaviour of Nx,y-i{X) is determined by nonzero singularities of o'x,y{t) , which are 
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much more difficult to study. According to Lemma [4.1.3[ the function Nx^y-i{X) 
can also be written as a convolution. The main technical problem is that the 
function pi has a jump, and therefore the straightforward integration by parts 
does not yield any new results. However, as we shall see in the next subsection, 
this jump disappears when we square N^^y-i^X) and integrate over y . Note also 
that |pi| is estimated by a smooth function whose Fourier transform has a com- 
pact support. This observation allows one to simplify and refine the well known 
Fourier Tauberian Theorems (see 18211. 



4.2. Upper bounds for Nx,y-i{X) . Since the eigenfunction are orthogonal 
in L'^{M), Lemma [4. 1.31 implies that 

(4.2.1) / \N,,y,,{X)\'dy = [ \ [ p,{\-ix)dNxM\'dy 
Jm Jm J 

pi{X-fi)pi{X-T)dN^^y{fi)dN^^y{T)dy = / \pi{X- dN^^M 

M J J J 

(an operator interpretation of this equality has been given in the subsection 12. 7p . 
The function |pi(t)P is continuous and infinitely differentiable outside the origin. 
Integrating by parts and changing variables in the right hand side, we obtain 



M 



N^,y;i{X)\^ dy = / |pi(A - p)pdA^^,^(p) = / pi{p) N^^^{X - p) dp 



where Pi(t) := ^ |pi(r)| is a rapidly decreasing odd function. Since pi is odd, 

/•OO 

Pi(p) N^,x{X - p)dp = J pi(p) {N^,x{X - p) - N^,x{X + p))dp. 
Substituting fll.2.2p in the right hand side and estimating 

i(A-p)"-(A+p)"i ^ c{\p\x^~' + \pn, \x±pr' ^ c(A"-i + ipr-i), 

we see that 

(4.2.2) / \N,^y,,{X)\^dy ^ CiX^-' + l), VA > , 

where the constant C depends only on the dimension, the remainder term in the 
Weyl formula fll.2.2p and the auxiliary function p. 

4.3. Upper bounds for A^j. y.o(A) . Since ^A^x,y;o(A) = / p{X — p) dNx^y{p) and 
p(A - p) dN^^yip) = p * iV;^(A) = {pa^^yY (A) , 



Condition 14.1.11 implies that the asymptotic behaviour of iVj. y.o(A) for large A 
is determined by the singularities of ax,y on supp p . For second order differ- 
ential operators, it is slightly more convenient to deal with the cosine Fourier 
transform ex,y(t) := f cos(tA) dNx,y{X) . The distribution ex,y coincides with the 
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fundamental solution of the wave equation. For all sufficiently small t, it admits 
the following Hadamard representation 

n+l 

(4.3.1) e.,y{t) = l^lE^^(^'^) r(j-Si + i) 

where Uj{x,y) are some smooth functions (see [Baj or |Berl Proposition 27]). Let 
us assume that 6 in Condition 14.1.11 is small enough so that fl4.3.ip holds on the 
interval (—6,6) (and, consequently, on suppp). 

Denote '$x,y{^) '■= Nx^y{X) — A^^ y(— A) . Changing variables and taking into 
account that the function p is even, we see that 

p * ^^,yiX) = N^,y;oW - iVa-,j/;o(-A) . 

Since Nx,y{\) = for A < and p is a rapidly decreasing function, Nx,y-o{—\) 
vanishes faster than any power of A as A ^ +oo. Thus we have 

(4.3.2) iVx,,;o(A) = P*^.,yiX) + o(A-"), VmGM+. 
Integrating by parts, we obtain 

(4.3.3) e,^y{t) = t j sin(tA) Ar^,j,(A) dA = -{2i)-Hi> .,^y{t) . 
In view of fHXT]) and (HX^ . 

(4.3.4) ^.,,(t) = -2zsignt^«,(x,y)i— -^(^ mod 
on supp p . We have 

(4.3.5) / e**^signt^ , "'^ , dt = 2i / sin(tA) ^ , , dt 

yiti>d,,, r(p+2) id,,, r(p+2) 



2 c^a-jj/ \ 
A 



J ^p(^dj-yX) 



(see, for example, [GRl formula 3.771(7)] or |GSl Chapter II, Section 2.5]). There- 
fore, taking the inverse Fourier transform of each term in the right-hand side of 
(14.3.41) . we obtain an asymptotic series 

oo 

(4.3.6) J2 ^) ^-'^^ "^f (^-.^^ ^) 

j=0 

with some bounded functions Wj{x, y). Let ^o-x,y{.^) be an arbitrary function such 

■ n 

that -^o;x,y{\) ~ 'LoWj{x,y)dx,y^ X'^ ^ Jii-j {dx,y X) as A ^ +oo uniformly 
with respect to x,y E M. By the above, 

(4.3.7) p * <ilx,yW - p * ^0;.,yW = 0(A-°°) . 
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Recall that the Bessel function Jpir) is estimated by C in a neighbourhood 
of the origin and does not exceed C for large r (see |GRt 8.440 and 8.451(1)]). 
Therefore 



(4.3.8) 



\dPX^Jp{d^^yX)\ ^ 



Cd: ' 



X > d~' 



x,y 5 



The inequalities fl4.3.8p imply that the terms in the expansion (14.3.61) are bounded 

■ n+l 

by Cj dx,y ^ X^ ^ for A ^ d~y and are estimated by Cj A""^-' for A G {e, d~y) , 
where e is an arbitrary positive constant. Thus, if A > e then 



^ Cj ( dx,y 



■I 2 



n — 1 

1 A— 



Wj{x,y)dly^ X"^ ^ Jn_j{da;,yX) 

for all j = 0, 1, 2, . . . and, consequently, 

— — n / ^ — ^ \ n — 3 

(4.3.9) ^0;x,yW-Woix,y)dxJ X^ Jn{d,^yX) ^ C[dx,'y + ij A" . 

n 1 

By direct calculation, wo{x,y) = 2~2 7r-2Uo{x,y) . We have 

uo{x,y) = {det{gjk{x,y)})~^^^ , 

where {gjk{x,y)} is the metric tensor in geodesic normal coordinates with origin 
X ([JB formula (3.1.1)]), so that Wo{x,y) = (27r)-t + 0(dj^^) (see [Roi p. 101]). 
From this estimate and fl4.3.9p it follows that 



M 



v^0;.,,(A)-(27r)-trf-|At Jn(4,,A) dy^QX^-' [d^-^ + l) dy 



AI 



for all X > e, where e is an arbitrary positive number and Ce is a constant 
depending only on e, the auxiliary function p and the geometry of M. Finally, 
since the function p is rapidly decreasing, the above inequality, (14.3.21) and (14.3.71) 
imply that 



(4.3.10) 



M 



Nx,y;oiX) - (271) 2 d:,J X"^ Jn{d^^yX) 



dy 



^ C, X""-^ [ (^,^7 + 1) dy , VA ^ £ > 

J M 



where is another constant depending only on e, p and the geometry of M. 
Remark 4.3.11. The above proof is a slight modification of arguments in |Baj . 
4.4. Proof of Theorem I2.1.2L The theorem follows from 04.3. lOp . K2:2^ and 

(03). 
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4.5. Proof of Theorem 12.3.11 By f l4.3.8p . there exists a constant Ci such that 



dx,y /i2 Jn (d^^yfl) 



^ Cl da;,' 



n + 1 



If C > V2C1 then the above inequahty and (12.1.31) imply that 



2 Cm ^ 2 



Nx,yi\) - C?x,y /i2 (4,2/ /i) 



dy 



N,^y{X)\'-2CU^r') 



4.6. Proof of Corollary 12.3.31 If y ^ x then the inequalities fl2.3.4p imply 
that y e Clk^N^^^i'^k, {.2C)~^Hk) for all sufficiently large N . By Theorem 12.3. ![ 
the measure of this intersection is zero. 

If /i^^ < 00 then, according to Theorem 12.3.11 the sum of measures of the 
sets f2a;(rfc, t/Ufc) is finite for all t > 0. By the Borel-Cantelli lemma, in this case 
almost every point y & M belongs only to finitely many sets Qxi^kitfJ^k) ■ This 

implies that limsup^^^o^ (^f^k^ ^x,yiTk)j < t almost everywhere. Letting t — > , 

we obtain (12.3.51) . 



4.7. Proof of Theorem Emu Let d 



M 



diam M . We have 



A/ 



x,y 



d,J A5 Jn (d^y A) dy ^ CX Ji {d,,y A) dy 

M ' 



^ CX I r''-^ jl{rX)dr 



CX 



XcIm 



r""' Ji(r)dr. 



Even if spherical coordinates centred at x do not exist globally on M, it is still 
possible to use the pull-back of the metric under the exponential map to change 
the integration variable to r. Since the Bessel function Js^ is bounded by C r 2 

in a neighbourhood of the origin and does not exceed C for large r, the 

right hand side is estimated by C A^^^ ^(n + + 

Thus we obtain 



AdM^^_2^^) if AciM > 1- 



d 



M 



x,y 



wo{x,y) dxj A2 Jn {d.^^yX) 



dy ^ 



C^{X'-- + d^i'), x^l 



CIlnAI 



x= 1. 



with = C\x — 1\ ^. This estimate together with (12.1.31) (multiplied by d^,^ on 
both sides) imply (12.4.21) by triangle inequality. 
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4.8. Proof of Corollary [2X3l Let ^^(A) := C^m (1 + A^-^) if x ^ 1 and 
g,{X) :=CM(l + |lnA|). By (1M2D, 

}Cxix,y)u{y)dy\ ^ g ^(X) \\u\\ 

M 



for all u e L°°{M) and 



M 



lCx{x,y)v{y) dy 



M JM \JM 



dx^ i }Cxix,y)dx)\viy)\dy^g^i\)\\v\\ 



for all V G L^{M) . These estimates imply the required result for p = oo and p = 1 
respectively. For other values of p the result follows from the Riesz interpolation 
theorem (see |So[ Theorem 0.1.13]). 

4.9. Proof of Theorem I2.2.1L If x > 1 then, in view of ([M2D, 

POO P 

/ / d:jN,,yi\)\'dydui\) < oo. 

Jo JM 

Therefore, by Fubini's theorem, the integral \Nx,y {X)\^ di'{X) is finite for al- 
most all y G M. 



5. Proofs of the lower bounds 

5.1. Proof of Theorem 12.5.11 Let the points x,y E M he not conjugate along 
any shortest geodesic segment joining them. Then, as follows from [Mill Theorem 
16.2], the number of shortest geodesic segments joining x and y is finite, and there 
exists e > such that {dx,y—e, dx,y+e)nCx,y = d^^y. Set ip{X) = exp(zA dx,y) p(A), 
where p satisfies Condition I4.1.1[ Then ip{dx,y) = 1, and one could choose p in 
such a way that supp ip C {dx,y — e, dx,y + e). By |LSVl Theorem 4.2] (see also 
|JPt Proposition 3.3.6]) we have: 
(5.1.1) 

(^^ax^yYip) = j i,{p-T)dNx,y{r) = C7e-'('^-^+*)/i^+0(/i^), 

where is the phase shift depending only on dimension of M (note that the 
Morse index of any shortest geodesic is zero). Since 

j W{p~T)\\Nx,y{T)\dr ^ I j i^\p~T)Nx,y{T)dT\ =1 j i^{p-T)dNx,y{T)\, 

(15.1.11) implies that 

j p'^\ilj\p-T)\\Nx,y{T)\dr ^ C/iV+'? + o(/i^+'^) , Vg>0. 
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The rest of the proof is similar to that of [Sari Lemma 5.1] Integrating the 
above inequahty and taking into account the estimate 



(5.1.2) / /iV+'^d/i » A^+^ 

Jo 

we obtain 

(5.1.3) y"^(A,r)|iV,,,(r)|dr ^ C A^+'' + O (a^+'') , 

A 

where ?/'(A,r) := J^^ jjf^ W{l^ — T)|d/i. Since i})' is a rapidly decreasing function, 
we have 



^(A,r) ^ C / (r^+ l/i-rH |V''(/i-r)|d/i ^ C (r^ + 1) , Vr G 



and '?/'(A, r) ^ (A |A — t|)~-' for all r ^ A and j = 1, 2, . . . These estimates 
imply, respectively, that 

(5.1.4) I i,{\T)\N^^y{T)\AT^ C f r'^|iV,,,(r)|dr 
Jo JO 

and 

(5.1.5) / V^(A,r)|iV,,,(r)|dr ^ C / v5(A,r)r"dr = 0(1). 
Putting together f l5.1.3l) -( r5.1.5p . we see that 

A-^-^ /'r«|iV..Jr)|dr » A"^""^ T r"^ |iV..,(r)| dr » 1. 



Jo 

Now (12.5.21) with p > 1 follows from Jensen's inequality. □ 

Remark 5.1.6. The inequality ( ]5.1.ip holds for Nx,x[}^ under the assumption that 
X is not conjugate to itself along any shortest geodesic loop. Therefore, similar 
lower bounds can be proved for the oscillatory error term R°^'^[X) introduced in 
\^V\ section 1.2]. Note that in dimension two this term coincides with the usual 
pointwise remainder in Weyl's law. 

5.2. Proof of Corollary 12.5.31 Condition (12.5.41) implies that there is a se- 
quence of points yUfc — s> +00 and a positive constant C such that 

Let Xfc(/^) be the characteristic functions of the intervals (0, /ifc] , A; = 1, 2, . . . By 
the above, the functions gk{f^) '■= lh:'^~^Xk{,lj) f^'^ {f{t^))~^ are uniformly bounded 
by a constant. Obviously, gk{f^) —^0 as — oo for each fixed /i. If the integral 
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/ /(/i) \Nx^y{fi)\^ dfi were finite then, by the Lebesgue dominated convergence 
theorem, we would have 



/i-^-^ / ^''|iv,,,(/x)rd/i= / <7fc(/i)/(/i)|iv.,,(/i)rd/i ^0, k 



oo . 



However, this contradicts to (12.5.21) . 

The estimates (12.5.51) and (I2.5.6P are obtained by taking /(/i) := and 
/(/^) := fJ-k^ i.xM - Xk-M) respectively. 

6. Examples 

6.1. Circle and 2-tori. In this subsection we justify Examples 13.2.21 and 13.2.41 

Let n = 1, 2 and let T be a lattice in M". Consider the flat torus T" = M"/r. If 
n = 1 one can assume that F = 27rZ, and thus = We have 

N.,yW = vor(n E/'^^^''"''^^^^"(2</A) 

where xb" is the characteristic function of the unit ball, Vol(r) is the volume of 
the torus and T* = e : (^, r/) G Z, V?] G T} is the dual lattice (see pi2| 
p. 29]. We are shghtly abusing notation here, since we added to the spectral 
function the constant term corresponding to ^ = 0. Note that N^^y = N^-yfi , 
hence we can assume without loss of generality that ?/ = . 

Let w{^) G C^(]R") be a nonnegative function depending only on |^|, such 
that / w(^)d^ = 1. Denote ^£(0 := e""' w{e~^^) for e > 0. By the Poisson 
summation formula, 

Voi(r) 



J2 f e2'^^<^+'''«>(ti;,*Xi?")(27re/A)de 
Tier -^K" 



A 



n 



^w{e\{x + T])) Xi?"(A(x + ?7)) . 



(27r, 

In what follows we set e = (AT)^^. Applying the asymptotic formula (11.3.31) for 

Xb" , one gets 



A^ 



2 A("-l)/2 sin (^A \x + r]\- / ;^(n-3)/2 



(27r)" ^ ^ (27r)("+i)/2 |x + r/|("+i)/2 ^ |x + r/|("+3)/2 

Using the formulae above, let us show that N^fi^X) G B'^ and 

2 sin ^A \x + ri\ — iZL^HzL^ 
(6.L1) iV,,o(A) ~ E (27r)(-+i)/2|a; + r/|(-+i)/2 " 
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The proof of flG.l.ip follows closely |B12[ section 3] and is split into four lemmas. 
Lemma 6.1.2. 



T 



dA 



A 



n-l 



0(1). 



Proof. Set 

/^(A) = XB"(27r^/A) - ^(AT)-! * XB"(2vr^/A) = Xxe'^C^tt^ - Wt-i * Xxb^C^t^O^ 

where A-B" is a ball of radius A. Since by definition J WT-^{C,)d^ = 1 
Wr-i £ C*^, it follows that /^(A) is supported on an interval of size 0(T^^) 
centred at A = 27r|^|. By a result of [Col] , we have 



#{e G r, lei < R} 



+ 0{R 



n-2+- 



yn = 1,2, 



Voi(r) 

Therefore, counting lattice points in an annulus of exterior radius |^| + 0(T^^] 
and interior radius |^| — 0(T^^) we get 



(6.1.3) e r*,supp/^(A) nsupp/5(A) ^ 0} = 0{T-'\^r-' + \^r'^^). 

This gives us an estimate on the number of non-vanishing cross-terms in the ex- 
pression under the integral in Lemma [6.1.21 Since all the terms in this expression 
are bounded by a constant, taking into account (I6.1.3P and the size of the support 
of /^(A), we obtain 



J2 e''^^^"'^^ (xB"(2vre/A) - ^(at)- * XB"(2vre/A)) 

27r|5|sCT 

Note that it is essential for the final equality that n < 3. 
Lemma 6.1.4. 



dA 
A^ 

-3+2/(n+l)\ 



0(1). 



□ 



T 



E 



IwiT ^ {x + ri))\ 



X\x + r7|("+3)/2 



dA 



0(1). 



Proof. Indeed, 



\w{T ^ {x + ri))\ 



X\x + r7|("+3)/2 



\w{T ^{x + ri))\ 



0(1). 
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Here the first inequality holds since A "^dX ^ 1, the second inequality is true 

(0 d^l ^ / ""^(0 = 1 5 cind in the last step we use 



because |w(a;)| = | / i 
that n < 3. 



□ 



In the lemmas below all the summations are also taken over elements rj eT. 
Lemma 6.1.5. 



lim lim — 



^ w{T ^{X + T])) . 



H>Q 



(n — 1 Ivr 



dA =0 



The proof of this lemma is a modification of the proof of |B12[ Lemma 3.3]. 
Proof. One can check that 



sin ( A |x + r^l 



in — 1)t{ 



sin I A |x + ,^1 



in — 1)t[ 



dA 



^ C min{T, I |a; + r^l — |a; + ^1 1 . 

Expanding the squared sum in the integral below and taking into account that 
|t&(a;)| < C(l + |a;|)~^" as w is rapidly decreasing, we obtain 



(6.L6) 



1 
T 



T 



E 

\r)\>Q 



w(T ^{x + r])) 



\x + ri\ 



sin I A |x + r^l — 



(n — 1)tt 



dA 



|-'?I>Q 



n-2+- 



|r7|('^+i)(l + T-i|r/|) 



4n 



c 



+ tH e 



M>Q k=\T-^\r,\'- 
+ 



\v 


1 ^ 

n + 1 


7] 


n-2+^ 

n + 1 


k 






1 + T-1 





n — 3 



(j + k) 



n — 3 
~2~ 



In the right hand side of f l6.L6l) . the first sum bounds the contribution of the 
terms corresponding to pairs ?7,^ such that ||a; + ?7| — + ^ . The second 
sum estimates the contribution of the terms such that T^^ < \ \x+r]\ — \x+^\ \ ^ 1. 

-. , 2 

Here we consider each sub interval of length |x + ?7| "+i separately. The last 
sum takes care of pairs 7], ^ such that k < | |a; + ?7| — |x + ^| | ^ k + 1, with /c ^ 1. 
For T > Q, the first sum in fl6.L6p is bounded by 



E 

\v\>Q 



T 


-1 


\v 


n—l _j_ 


V 


n 






r]\ 


"+1(1 + T- 


1 





C {Q-' + T-'Q-') , 

C {Q~^/^ + T~^\\nT -\nQ\) 



n 



n = 2. 
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In order to estimate the second sum in (16. 1.61) . we first note that 



M 



-J 



fc=[T-i|r?|^ ~ 



1^1 



< C\r]\^~^ InT. 



for all sufficiently large T. Estimating the sum over t] by an integral we get 



\v\>Q 



\v 


1 

11+ 1 


V 


n-2-\ — ^ 


k 


V 


n+1 ^ 


1 + T-i 





fc=[T-i|r;|""" 



< 



ClnT 
T 



+CO 3 



1 + r-ir 



CqT-MuT, n = l, 



^ n— 3 n— 3 

To find a bound for the third sum in (16.1.61) . we use the inequality 2a^~b^~ ^ 
_j_ fjn-3 gj^^^ oucc again estimate the sum by an integral, 



J 2 (j + fc) 2 



+ r-ij7"(i + T-i(j + fc)) 



2n 



< 



c 



hoo „n— 3 



Hence we have 

hoo 



J 2 (j + /C) 2 



+ c 



Q+k 



oo ^n-3jj^j, 

(1 + T~ir)2" 



- — 

— 1 j^Q 

This completes the proof of the lemma. 



CgT-MnT, n = l, 
CnT-Mn'T, n = 2, 



□ 



Note that the proof of Lemma 16.1.51 uses the fact that the limit with respect 
to T is taken first. 



Lemma 6.1.7 



IvKQ 



1 - w(T-\x + ri)) . / , , (n- l)n 
sm A |x + ?7| — 



|a; _)_ ^|(n+l)/2 



4 



dA = o(T). 



Proof. This estimate holds because the sum involves a finite number of terms and 
limT_^+oo ""^(^"^a^) = 1 for all x. □ 



By Lemma [6.1.21 we get 



lim — , 

T^+ooT 



2 dA 

\n-l 



0. 
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Lemmas I6.1.4[ 16.1.51 and 16.1.71 imply 
1 



lim lim 



SO that 



((AT)- 



xfi 



2 sin Ix + ?7| — lliJ^ j 
(27r)("+i)/2 \x + r7|("+i)/2 



dA = 0, 



lim lim — 



2 sin ^A |a; + ?7| — 

(27r)("+i)/2 \x + r/|("+i)/2 



(n— l)7r 
4 



dA 



0. 



This implies (16.1.11) and therefore proves Examples 13.2.21 and I3.2.4[ 

Remark 6.1.8. As was indicated in Remark 13.2.61 this approach can not work for 
n ^ 3. Indeed, in higher dimensions we can no longer bound the sums under the 
integral in Lemmas I6.1.2[ 16.1.41 and 16.1.51 by taking the absolute value of each 
term. A more delicate analysis is required in this case (see |BB] for some related 
results). 

We also note that there is a simpler and more direct proof of Example 13.2.21 
based on the identity \GR\ formula 1.422(4)] 



(6.1.9) 

It works in dimension n 
reader. 



1 



sm 



2/S-' 



E 



{s + 2nky 



2^ fc=- 

1 only, and we leave the details to the interested 



6.2. Spheres. In this subsection we are going to justify Examples ll.3.4l and l3.2.7[ 
The spectral function on depends only on the distance or, equivalently, on the 
angle s between the points x and y. For any x G S", the eigenfunctions orthogonal 
to the Legendre polynomials Pm{n,t), with t = coss and m ^ 0, vanish at x. 
Hence, Pm{n,t) are the only eigenfunctions that matter to compute the spectral 
function. They satisfy the differential equation 

((1 - t^)d^ - ntdt)P^{n, t) = -m{m + n - l)P™(n, t) , 

which is equivalent to the eigenvalue problem A/ = m(m + n — 1)/ on with 
/ depending only on the parameter t G [—1, 1]. Note also that 



PlXnM^-e^'^dt 



1 



Pm{n^) = 1, 



where 



r(H±i) 



and N{n, m 



Dn-i N{n,m) 

(2m + n — l)r(m + n — 1] 



r(m+ l)r(n) 
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(see |Mu[ Lemma 10]). We assume that x and y are non-conjugate points, so 
that < s < TT. The rescaled spectral function on S", for n ^ 2, is given by 

N{n, m)Pm{n, cos s) 



D. 



n 



0^m(r?i+n— 

Here we are abusing notation shghtly, since according to (11. 2. II) the inequahties 
in the summation hmits should be strict. However, because of the rescaling this 
makes no difference m. B^. 

Consider the following generating function ( |Muj . Lemma 17): 

+00 2 

(6.2.1) J]iV(n,m)P^(n,t)z™ = 

[1 + z'^ — 2zt) 2 

Zk=o^kZ^ = /(^) in 

Therefore (16.2.11) implies 



m=0 

It is easy to see that Yl'k=o'^kZ^ = fi^) implies (Z]i=o%) = 



(6.2.2) ^A/',(v/m(m + n- 1))^™ 



Dn{l + - 2zt)^ 



m=0 

Lemma 6.2.3. For m E N, the spectral function on satisfies the following 
asymptotic formula 



Ns{\/m{m + n — \)) 



n — 1 

2cos(|)m~ f fn ^ \ (n + 1) \ n-s. 



- cos - + m s TT +0(m2), m oo . 

(27rsins) 2 \\l / 4 / 

Proof. Following the approach of [CoHt Chapter VII, Sections 6.6 and 6.7], we 
prove Lemma [6.2.31 using the Darboux method applied to the generating function 
(16.2.21) . The Darboux method relies on the following fact: the coefficients of 
the Taylor expansion at the origin of the function f{z) = 'Y^^=q cikz'^-, holomorphic 
in the open disc D, decay as 0{k~^) if /(e*^) G C""(M). The first step to obtain 
the asymptotic formula is to approximate the generating function (16.2.21) . taking 
into account the singularities of the highest order. 
Near the singular point x = e^**, we can write 

n+l 

1 {{e^'' - e^'') + {z - e^'')) ^ 



z^-2zt)— {z-e^'') — 



e^-^ ^ (-^\ (z-e^ 

(2sins)^(z-e±-)^ V k )\±2i^ms 



k=0 

where 

'a\ I 1, k = 
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For p > we have: 
(6.2.4) 



(2sins)'^^V k jy (2^sins)*^ (-2isins)'= 

^ ^ e 

Expanding the powers of 2; — e^*^ by Taylor's formula at z = 0, we see that 
Taylor's expansion at z = of the difference (16.2.41) is given by 

2 Y^Y^/-^\ A-^ yos((^ + /-fc). + (/-^-| )7r) , 
(2sin.)^^t^V k )[ I ) (2sin.)'= 



2r . 



Using (16.2.21) and comparing the coefficients in front of the same powers of z, we 
deduce that 



DnNs ( y/m(m + n-l)) 

_ 2 >A f-^\ ( cos ((^ + m - fc)s + (m - ^ - |)7r) (k-^ 

n + l 



(2sin.)^ V ^ ; V (2sins)'= 



m 



k=0 

2 _ ^ /^-^^ cos ((^ + m - fc)g + (m - ^ - |)7r) ^A; - ^ 



(2sins)'^^V k jy (2sins)'= V " 1 

+ 0(m~i). 

as m ^ 00. Here the error estimate follows from (16.2.41) in view of the remark in 
the beginning of the proof. Taking into account the estimate 

1 (m + 1) 

we see that the main term of the expansion above corresponds to k = 0. Hence 
we obtain 

D^Ns ( ^/m{m + n-l)) 

4cos(f)r(m+^) //n \ (n + l) \ ^, 
- ^ ^ cos - + m s-^ ^TT +0(m 2 



(2sins)^r(m + l)r( 
This completes the proof of the lemma. □ 
Note that on S" there are simple expressions for Ns{fi) at conjugate points: 

{n + 2m) {n + m - 1)\ 27r"/2 



A'„(Vm(m+ «-l)) = = p^)..r(„^a)„ ™ + 0(m , 
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These formulae and Lemma 16.2.31 prove Example 11.3.41 

Let us now complete the proof of Example I3.2.7[ Lemma 16.2.31 implies that 

2 

dX 



2cos f //^ n-1. n\ {n + 1 

Ns{X) - , ^ cos MA - + 77 s - 



(27rsins)' 



2 ' 2 I 4 

T 



< c$:i=o(i„T), 



m 

m=l 



Here we have used that the interval 

^ym{m + n-l), ^/ {m + 1) {m + n)^ p| {A : [A - (n - 1) /2j 7^ m}, 

where the formula above does not agree with Lemma r6.2.3l is of size 0(m^^). This 
follows from the simple asymptotic formula ^Jmim + n — 1) = m+^^+0(m^^). 
Therefore, Ns{X) is 5^-equivalent to 

2cos(f) // n-1 n\ {n + 1) 
— cos MA —J + -]s TT 



(27rsins)^ VV 2 ^ 2^ 4 

which can be rewritten as 

1 1 n — 1, n\ (n — 1 

(6.2.5) ^^^-^.^^^^^ sin Ma - 



(27r)^sinV(s)sin(f)"- VV"' 2 ^ ' 2/ " 4 
Lemma 6.2.6. The following expansions hold in for < s < vr.- 



-TT 



sin((LAJ+i).) 



+00 



2sin(f) ^ |s + 27rA; 

2sm(^) s + 27rA; 

^ fc=— 00 

2sin(|) 5: ^T^^°<^l^ + 2-^l)' 

^ fc=— 00 

COs(LA+ij5) ^ i^-lf+Hik) 

— ~ > n — rr COS A s + 27rA; ) . 

2 sin I ^ \s + 2Tik\ ^ ' 

^ fc=— 00 

HereH{x) is the "reversed" Heaviside function: H{x) = ifx ^ 0, andH{x) = 1 
ifx<0. 

Proof. Each of these expansions can be obtained in a similar way as fl6.1.ip in 
the case n = 1. Alternatively, they could be proved directly using the identity 
f l6X9D . □ 
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Applying Lemma [6.2.61 to the formula fl6.2.5p . we deduce that 

oo sin(A/(7,)-^-c.(7,)f) 



TT (27rsins) V ^ l{-fj) 

where the sum runs over all geodesic segments 7^ starting at x and ending at 
y. Here l^jj) is the length of jj and uj{'jj) is the Morse index of 7^. Counting 
the points conjugate to x on the geodesic 7 of length |s + 27ik\, one gets uj{'y) = 
{n — 1) (2|A;| — H{k)). This completes the proof of Example 13.2.71 

7. The (-function: proofs 

7.1. Auxiliary functions. Throughout the section we denote 

• t:=Yiez, s:=lm2; and 

• (s) := (1 + |s|2)i/2. 

Let p be a real- valued even function satisfying Condition 14.1.11 with a suffi- 
ciently small 5 . Let us fix an arbitrary c > such that 2 c < Ai and c < d~^y 
for all x,y E M , and denote 



(7.1.1) Z^^y,^{z) := z J X'^-' N^,y;jW dA , j= 0, 1, 
where N^^yj^X) are the functions defined by (14.1.20 . Then 

and, due to the finite speed of propagation, \Zx^y{z) — Zx,y-i{z)\ ^ Ct whenever 
5 < d^,y ■ 

By Lemma [4.1.3[ we have 

(7.1.2) Z^^y.iiz) := [ hiz, fx) dN^M , 



where hi{z,fi) := z X ^ ^ pi{X — fi) dX . Substituting zX ^ ^ = — ^ (A ^) , 
integrating by parts and then changing variables, we obtain 

POO 

(7.1.3) hi{z,n) = + pi{c- n) - (A + ;u)~XA) dA, ^ii>c. 

J C—fJ, 

Since one can differentiate under the integral sign, hi {z, /i) is an entire function 
of the variable z E C smoothly depending on the parameter p G (2c, +00) , such 
that 

d"^ T(z-\-m) 

(7.1.4) 1^ h{z, p) = {-ir -^YW^ hi{z + m, /i) . 

Lemma 7.1.5. For all p G (2c, +00), we have \hi{z,p)\ ^ Ct and 
(7.1.6) \hi{z,p)\ ^ a,^/i'*'"(s)"(l + |ln(s) -lnpr^ Vx,r>0. 
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Proof. The identities (17. 1.31) and / p(A) dA = 1 imply that 



hi{z, /i) = c ^pi(c - /i) + p" 



p(A)dA+ / (p-^ - (/i + A)-") p(A) dA 



+ 



a 

2 



(p-^ - (p + A)-^) p(A) dA + / (p^^ - (p + A)-^) p(A) dA 



Since p is rapidly decreasing, the first four terms in the right hand side are 
bounded by Ct^mfJ'~"^ for all m > 0. Therefore we only need to estimate 

/j^((p + A)--p-^)p(A)dA. 

By Taylor's formula, 



m—l 



(7.1.7) (p + A)-^ - p-^ = J2 (-1)' ^ffw,? -^V'"' 



fc=i 



T(z + m) 



dt 



m = 1,2, 



(m-l)!r(2;) 7o 

(unless —2 is a nonnegative integer, in which case the terms with negative expo- 
nents —z — k are absent). Substituting fl7.1.7p into the integral, we see that 



(7.1.8) / ((p + A)-^-p-^)p(A)dA 



m—l 



(-i)^r(^ + fc) 



+ (-i: 



A;!r(2) 

Viz + m) 



A'=p(A)dA 



(m - l)!r(z) 



(1 - r)"'-^ (p + rA)-"-™ A™ p(A) dr dA . 







Since p is rapidly decreasing. Condition 14.1.11 implies that the sum in the right 
hand side of 07.1.81) is estimated by CmP"*"™" {s)^~^ for all m ^ 1. On the other 
hand, for all j, m ^ , 



(7.1.9) 



|(p + rA)-^-'"AV(A)|dA ^ a,„p 



-t—m 



Vr G [0, 1] 



Therefore the last term in (17.1.81) does not exceed Ct^m P * (s)™". Thus we have 



(7.1.10) 



(p-^-(p + A)-^) p(A)dA 



t im 



Vm = 0, 1, 



where bg^^ := p ^{s) (the estimate with m = is a particular case of (17.1.91) ). 
Since (1 + | Infe^^^l)"'' ^ Cr,^ min{6s,^''^, b7/i^} whenever x > and 

p- (1 + I ln(.) - Inpl)-'- = 6-. (1 + I ln6,,,|)-^ , 
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interpolating between the estimates fl7.1.10p . we arrive at (17. 1.61) . □ 

7.2. Estimates for Zxy.i{z). Lemma [7.1.51 implies the following corollaries. 

Corollary 7.2.1. Zj.^y.i{z) is an entire function of z for each {x,y) G M x M, 
such that 

(7.2.2) \Z,^y.,{t + ts)\ ^ + yt^n. 

Proof. In view of (17.1.4^ and f l7.1.6p . one can differentiate under the integral sign 
in the definition of Zx^y.i{z). Therefore the function Zx^y-i{z) is analytic. 

Let Nx,y{\]ai,a2) ■= J2xj<x\^iVji^) + V^i(2/))P where ai,a2 are complex 
constants such that + |a2p = 2. The function N^^y{X; Oi, 02) is nondecreasing 
and, in view of the Weyl formulae (ll.2.2p . is estimated by C A". We have 

2 ReiV,,j,(A) = N,,y{\; 1, 1) - iV,,,(A) - Ny^y{\) , 

2t ImN^^yiX) = N,^y{X; t, 1) - iV,,,(A) - Ny^y{X) . 

Therefore it is sufficient to prove that the estimate (17.2.21) holds for the integral 
I2C ^) dG'(A), where G(A) is a nondecreasing function bounded by C A" . 
U t < n then (I7.1.6P with x = n — t and r = 2 implies that 

/ hi{z,X)dG{\) ^ Ct{s)''~' A-"(l + |ln(s) -lnAI)-MG(A) 
J2C J2C 

POO 

^ Ct (s)"-* / A-"-^ (1 + I ln(s) - In A|)-2 G(A) dA 

POO 

/ (1 + |ln(s) -InAD-'A'MA ^ Q . 
J2C 



^ Ct{s) 



li t > n then the required estimate is obtained in a similar way, with the use of 
the inequality \hi{z,\)\ ^ A~* . □ 

Remark 7.2.3. Note that (17X2]) holds for x = y. 

Corollary 7.2.4. 

(7.2.5) / \Z,,y,i{t + ts)\'dy ^ Ct{{sr-''+l) , \ft - , \/x e M . 
Proof. In the same way as in (14.2. ip . we obtain 

JjZ,,y,,iz)\^ dy = J |/ii(z,A)|MiV,,,(A). 
If t < n/2 then (17.1.60 with x = | — t and r = 2 implies that 

/oo 
A-"(l + |ln(s)-lnA|)-MiV,,,(A) 

POO 

^ Ct{sy-^' / A-"-i(l + |ln(s)-lnA|)-2iV,,,(A)dA. 
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Using the Weyl formula and estimating integrals as in the proof of the previous 
corollary, we see that the right hand side is not greater than Ct . 

If t > I then the required inequality is obtained in the same way, with the use 
of the estimate A)p ^ A'^*. □ 

7.3. Estimates for Zx^y-o{z). \i x ^ y then Zx,y-o{z) is an entire function of 
z. Our goal is to estimate Z^^yfi^z) uniformly with respect to s and x,y E M . 
Further on 

• Ot{d^y) denotes a function which is estimated by Ct{d^y + 1) . 

Obviously, if |/(A)| ^ C(A)-™ for all A G R+ then J^X~^f{X)dX is an 
analytic function in the half-plane {z : t > 1 — m} and A~^ /(A) dA| ^ Ct 
for all t > 1 — m. In particular, if \E'2,._j,(A) = Nx,y{X) — Nx^y{—X) then, in view 
of (1133), we have 

/"OO /"OO 

Z.,y;o{z) = z / A-^-1 N^,y;oW dX = / A-^ (p * ^^J (A) dA + 0^(1) . 



Let us assume that 6 in Condition 14. 1.11 is small enough, so that the Hadamard 
representation (14.3. ip holds on supp p. Then the Fourier transform of the deriv- 
ative "^'^ y admits an asymptotic expansion of the form (14.3.11) on supp p . This 
implies that 



/oo 
P(A - /i) d/i = 0{X 



uniformly with respect to x,y E M , where ^E'q.^. ,^(p) is a function whose asymp- 
totic expansion for /i — ^ oo is obtained from (I4.3.6P by differentiating each term 
with respect to A and putting X = p . Thus we have 

/oo /"OO 
J X-^piX-p)%^,^yip)dpdX + Oi(l). 

For all m = 0, 1, 2 . . . , we have Jpi^) ^ C tP~^ for small values of r and 
-^Jpir) ~ r"2 cosr^a^^fcT"'' + r"2 sinr^6„,fer"'' , r ^ oo , 

k=0 k=0 

where ^ and bm,k are some real coefficients (see, for example, |GRt 8.440 and 
8.451(1)]). Therefore the asymptotic formula for "^'o-x^y implies that \'^o-x,yif^)\ ^ 
C (p)"~^ for all p e (c, d~y) and 



(7.3.2) = 5Z %,fc(a:,2/)4'ir'""«yP)'^"'"'exp(24,,p) 

j+k<N 

+ X] c^Jy^^"" (4,j/P)^ ^'"^ exp{-idx,yp) + RNix,y,p) 

j+k<N 
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for all II e (^xi' ^) N — 1,2, ... , where Uj^k , ^ ^-re bounded functions 
and 

\RN{x,y,ii)\ < Cdl-^ (4,./^)"^-^ , V// e Kj,oo) . 

rt— 1 

In particular, |^o:x',y(Ai)| ^ C'c^i.^" (c^rr.y/u) ^ ^ C fj, whenever fj, ^ d~y and, 
consequently, |^o;j:^y(/x)| ^ C {^)^^^ for all /i > c. 
The above estimates and the obvious inequalities 

(7.3.3) {X)-\^x) ^ 2(A - ^i) , (//)-i(A) ^ 2(A - ^i) 

imply that 

\x-'p{x-t^)%,^,M ^ a (Ar-*-^(A-/.ri|p(A -/.)!, 



and 



I I , . J. n- 

A"V(A-/i)i?iv(a;,y,/x) ^ C^c^L (4,yA) p(A - ^) (4 .^/x) ^ 



t-n+l 
'x,y 



{dx,yX) ■ 



-N-t 



for all X, II > d^ y. Since the function p is rapidly decreasing, integrating the 
above estimates over A and /i, we obtain 



(7.3.4) 

(7.3.5) 

and 

(7.3.6) 



X 'P{X- ll)%.^^y{fi)dfidX 



X-^p{X-ti)%.^,^y{ti)dtidX 



'a(C/ + i), t^n. 



CWnd 



t — n 



Ct (rf*-" 
Cllnd 



■x,y\ 1 



'x,y \ 1 



t = n 



poo poo 

/ ^ / ^ A^XA - p) Rn{x, y, ii) d/i dX 

*^ , y , y 



OO POO 



< C,d 



t-n-1 




^ ad 



-N-t 



-N-t 



dX 



{X-ii)'^^''\p{d-i{X-ii))\diidX 



for all N > ^-t. 



It remains to estimate the integrals 

poo poo ^ 

(^■3-7) _J_^ ^"V(A - p) dg+'-^ {d,,yp)'^-^-' e^'^^-'^dpdX 

) 

X-^p{d-l(X-p)) p'^-^-'e^^'^dpdX 



j2j-i-t-n-l-\-is 
"'x,y 
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generated by the main asymptotic terms in (]7.3.2p . Expanding fi~ ■'^^ by 
Taylor's formula at /i = A, we see that the right hand side of (17.3.71) coincides 
with 

/OO / POO \ 

^l^iz, A, /x) ^ p {d-^li\ - /i)) e±^^ d/ij dA + 0,(rfj^*-") , 

where 

I / 72— 1 ■ / \ 

^(2,A,/i) = ^ / — j A'^-^-'^-'"-^ (/i - A)'" 

m=0 \ ^ / 

with an arbitrary Z ^ — t — j — /c , and the remainder estimate follows from 
the inequalities 

and 

Note that the last inequality is a consequence of (17.3. 3p . 
Obviously, 

OO r-\ 



d-'y / / (/i - A)™ |p (rf-,^(A - p)) I d/xdA 

J 1 J — OO 

poo POO 

^d^'J / An/i™p(/i)|d/idA ^ Cp,^, VJ9GM, Vm = 0,l,2,..., 

Ji Jd-}(\-i) 

where (Im is the diameter of M . This estimate and (I7.3.3P imply that the 
integral (I7.3.7P is equal to 

/OO / POO \ 

V^(z,A,/x) pKJ(A-/i)) e^^'^d/ijdA 

V 2 ^ / A'^-^-'^-^-^e^'^dA 

m=0 V / ' ' ^1 



modulo Ot{dl^^^ ") , where the mth derivative of the Fourier transform. 

For each p G M , the integral A^^^ e^"^^ dA defines an analytic function on 
the half-plane {z : t > p + 1} , where it is bounded on each vertical line by a 
constant Ct ■ This function admits an analytic continuation to the whole com- 
plex plane, obtained by replacing e^^^ with (=fO"^ ^^^^ ^^^1 integrating by 
parts. This continuation coincides with the difference between the meromorphic 
continuations of the integrals \^~^ e^*'^ dA and A^^^ e^*'^ dA . According 
to [GEl 3.381] and [GRl 8.328(1)], 

POO 

/ A^-^e^^MA = ±ie^'^^^-'^'^T{p + l- z) 
Jo 
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and 



whenever |s| > 1. Replacing ^ with {m+p- 



■z 



t + is 

-1 



|p-t+l/2 



and integrating by parts, we see that 
Therefore 



dA 



l_d™_ \m.+p—z 



^ Cp_t whenever \s\ > 1 . 



AP- 



-2: _±iA 



dA 



p-t+l/2\ 



The above inequahty imphes that the integral in the left hand side of (17.3.81) 
is estimated by C't(i|{+*~" (1 + (s) 2"~*~^~'^) and, consequently. 



(7.3.9) 



1 

'x,y 



A-V(A-^)<^i-"(4,,/^)^ 



x,y 



dfi dA 

(l + (s)t-*-^-'=) 



Now, putting together (I7.3.4p -( 17.3.6p and (I7.3.9p . we obtain 



(7.3.10) 



C((.)t-*+|ln4,,|) 



7.4. Proof of Theorems 12.6.51 and I2.6.6L Theorem 12. 6. 51 is an immediate con- 
sequence of (I7.2.2P and (17.3. lOp . Since the function (c?^*"""^ + l) is bounded 



and is estimated by d^ y 
dZjSD and (17.3.101) . 



n-2t+e 



for small values of d^^y , Theorem 12.6.61 follows from 
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